Form-factors in the Baxter— Bazhanov— St roganov 
model I: Norms and matrix elements 



! G von Gehlent, N Iorgov*, S Pakuliak*, V Shadura* and 

Yu Tykhyy* 

. t Physikalisches Institut der Universitat Bonn, Nussallee 12, D-53115 Bonn, Germany 

O * Bogolyubov Institute for Theoretical Physics, Kiev 03680, Ukraine 

" Bogoliubov Laboratory of Theoretical Physics, Joint Institute for Nuclear Research, 
Dubna 141980, Moscow region, Russia 



CO 



(N ■ Institute of Theoretical and Experimental Physics, Moscow 117259, Russia 

E-mail: gehlen@th.physik.uni-bonn.de, iorgov@bitp.kiev.ua, 
CO' pakuliak@theor.jinr.ru, shadura@bitp.kiev.ua, tykhyy@bitp.kiev.ua 

d 

Abstract. We continue our investigation of the Zjv-Baxter-Bazhanov-Stroganov 
model using the method of separation of variables [T]. In this paper we calculate 
the norms and matrix elements of a local Z^r-spin operator between eigenvectors of 
\ the auxiliary problem. For the norm the multiple sums over the intermediate states 

' are performed explicitly. In the case N = 2 we solve the Baxter equation and obtain 

form-factors of the spin operator of the periodic Ising model on a finite lattice. 
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PACS numbers: 75.10Hk, 75.10Jm, 05.50+q, 02.301k 
1. Introduction 



The Baxter-Bazhanov-Stroganov model (BBS-model, also known as the r^-model) 
[21 [3l HI [5] has attracted considerable interest, because via functional relations it is 
related to the solvable ZAr-chiral Potts model (CPM). Solving these functional relations 
has been the main method to calculate the eigenvalues of the CPM [6J . However, also by 
itself, the BBS-model is an interesting lattice spin model with cyclic Z^r spin-variables, 
which is closely related to the six- vertex model at roots of unity [THUE]. The BBS model 
in its vertex formulation can be solved using the Functional Bethe ansatz or Separation 
of Variables (SoV) method [9J. For iV = 2 it is equivalent to a generalized free-fermion 
Ising model [ID] . 

Using the formulation of the BBS-model in terms of cyclic L-operators [3, HUH], in 
[I] we developed a general method for the construction of its transfer matrix eigenvectors. 
Our approach is an adaptation of the SoV method of [9] [12] to iV-state spin chain models. 
In order to find the eigenvectors of the transfer matrix we first obtain the eigenvectors 
of a certain auxiliary integrable system. Commuting integrals of this auxiliary system 
are generated by the off-diagonal elements of the monodromy matrix of the BBS model. 
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Then the eigenvectors of the auxiliary system serve as building blocks for the eigenvectors 
of the periodic BBS model. The multi-variable kernel which relates the respective 
eigenvectors can be presented as the product of single variable functions (SoV). Each 
of these functions satisfies the Baxter equation. In [Tj we showed that the existence 
of non-trivial solutions of these equations is equivalent to the well known functional 
relations [21 HI E] for the transfer matrix of BBS model. 

The goal of this paper is to use the explicit constructions of p] for calculating matrix 
elements or form factors of local operators of the BBS model. Since the eigenvectors 
are given in terms of multiple summations, the main concern of our calculations is to 
perform these sums explicitly in order to get factorized expressions for the form factors 
and norms. For N = 2 from the results of our first paper we calculate spin operator 
matrix elements of the generalized Ising model. 

We hope that the methods developed in [1] and in this paper will allow to find 
analytical formulas for matrix elements of the CPM. 

We would like to mention the paper [13] , where matrix elements of local operators 
for the Toda chain are calculated using the eigenvectors of [T2] . 

This paper is organized as follows: In Section 2 and 3 we recall the results of pQ 
for the right eigenfunctions of B n (X), adding the analogous results for the respective 
left eigenfunctions, pointing out the differences. Then we calculate the overlap of the 
right and left eigenfunctions of B n (X), i.e. the norms. Here the main result will be that 
the sum over intermediate quantum numbers can be performed explicitly. In Section 5 
we concentrate on the case N = 2 and calculate the overlap of the eigenvectors of the 
periodic transfer matrix t(X) = A n (X) + D n (X) for the homogeneous BBS model. At 
the end of Section 5 we also write a result for the matrix elements of the single particle 
operators u n which still has multiple summations over discrete internal variables. In 
Section 6, using explicit solutions of the Baxter equation for the case of the homogeneous 
Ising model, we announce the factorized expression for the matrix elements of the one 
site spin operator where all internal summations have been performed. The details of 
this calculation will be given in a forthcoming paper [16]. These will prove a result 
conjectured by A. Bugrij and O. Lisovyy for the matrix elements of the spin operator 
in the Ising model on a finite lattice [HI [15] . Finally, Section 7 summarizes our results 
and the Appendix explains our summation technique. 

2. L-operator formulation of the BBS-model 

We define the BBS-model as a quantum chain model in the L-operator vertex 
formulation (the relation to the formulation as a face model can be found in [HE]). 
To each site k of the quantum chain we associate a cyclic L-operator [7J HJ [TT] acting in 
a two-dimensional auxiliary space 
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At each site k there are ultra-local Weyl elements and obeying the commutation 
rules and normalization 

S ■ i 2iri IN N N i 

uyUfe = u k uj , VjV k = v k Vj , UjV k = uj J - k v k Uj , u) = e 7 , u fc = v fc = 1 . 

A is the spectral parameter. We have five parameters x fc , a k , b k , c k , d k per site. At 
each site k we define a iV-dimensional linear space (quantum space) V k with the basis 
\j)k, 7 £ Zjy, the dual space with the basis fc(7|, 7 £ 2W, and the natural pairing 
fc(7'|7)fc = 5y )7 . In Vfc the Weyl elements and act by the formulas: 

u fc |7)fc = w 7 |7)fc , v fe |7)fc = |7+ l) fe . 

In VI these Weyl elements act as follows: 

fc(7|u fc = fc (7k 7 , fc(7|v fc = fe (7 - 1| • 

The monodromy matrix for the quantum chain with n sites is defined by 

r .(A) = MA) *<*)■■• MA) ^J). (2) 

The transfer-matrix is defined to be the trace in the auxiliary space t n (A) = 
tr T n (X) = A n (X) + D n (\) . This quantum chain is integrable because the L-operators 
(PQ) are intertwined by the twisted 6-vertex i?-matrix at roots of unity. It leads to 
[t n (A), t n (/i)] = and so t n (A) is the generating function for the commuting set of 
non-local and non-hermitian Hamiltonians Ho, . . . , H n of the model: 

t n (A) = H + H l X +••• +H n _ 1 A"- 1 + H n A". 

From the intertwining relation, it follows that B n (X) is the generating function for 
another commuting set of operators hi, . . . , h„: 

[B n (X), B n (jj,)) = 0, B n (X) = h x A + h 2 A 2 + • • ■ + h n A" . 

Following Sklyanin [9], we shall first solve the eigenvalue problem for this last set of 
commuting operators, which turns out to be possible by a recursive procedure. After 
this, the periodic problem will be related to this auxiliary result by Baxter equations. 
After proper normalization by a Sklyanin's measure, the kernel of the Q-operator 
factorizes into a product of single variable functions (SoV). 

3. Eigenvectors of B n (X) 

3.1. New basis vectors, one-site states 

In the root of unity case an important role is played by the cyclic function w p (j) [TF] 
which depends on the Z^-variable 7 and on a point p = (x,y) on the Fermat curve 
x N + y N = 1. We define w p (j) by the difference equation 
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The Fermat curve restriction guarantees the cyclic property: u> p (7 + N) = Wpi^y). The 
function Wpipf) is a root of unity analog of the g-gamma function. 

It is convenient to change the bases in the spaces V k and V k . Instead of \^) k and 
fe(7|, 7 G Z N , we will use the vectors 

\^p)k = ^2^ P f(i - p)h)k, k(ip P \ = Y^ — T~Z — T~n fe ( 7 '' pe^iv- (4) 

The coordinates of the Fermat curve points p\ = (x k ,y k ) and = (x^,yf) are 
defined as follows. Let us fix some value of r k to satisfy = a k — b k . (We shall 
consider generic parameters such that r k 7^ 0.) Then 

L R r k R b k l ^ 

X A; — X k — X k — , Uk — Vk — > Uk — • 

a k a k tua k 

Observe that yf = ujy\ , while x\ = xf. The vectors (jlj) are eigenvectors of the upper 
off-diagonal matrix element Xu k 1 (a k — b k v k ) of the operator L k : 

Xu k 1 (a k - b k v k ) \if) p ) k = Xr k u~ p \ip p ) k , (5) 

k (il) p \Xu k 1 (a k - b k v k ) = \r k u~ p k (4> p \ . (6) 

The first equality has been proved in [1] , the second relation follows very similarly. The 
action of the operator is 

fc(Vv|vfc = k{ip P -i\ , VjtlV'p)* = \tp P +i)k ■ 

3.2. Right and left eigenvectors of B n (X) for general chain length n 

The eigenvectors of B n {X) are labelled by the vector p n = (p n ,o, ■ ■ ■ , p n ,n-i) £ (2W) n - 
Let us further define: 

Pn = Ylk=0 Pn,k] Pn = {Pn,l, Pn,n-l) £ (Z^)™" 1 

and p^ k denotes the vector p n in which p n>k is replaced by p Hjk ± 1, i.e. 

Pn k = (Pn,0, Pn,k ± 1, • • • , Pn,n-x), fc = 0, 1, . . . , 71 - 1. 

The following formulas give an iterative procedure to obtain the eigenvectors of 
B n (X) from eigenvectors of B n _i(X) and single site vectors defined by (jlj). 
The vector 

\^ Pn )= J2 Q R (pn-l,Pn\pn)\V Pn - 1 )®\lp Pn )n (7) 

Pn-ieCZiv)™- 1 

Pn 

, /.,(pn-pn-l)(pn-Pn,o) 

^ ^^(Pn-l, P>n) = ? " rr— p 7T X 

«V,o(Pn,0 - Pn-1,0 - ij^PnlPn - Pn ~ 1) 



X 



nr=i 2 re; ^j^ 1 - p ^ k) g u^h^: 



n^=i W „»-1> 1 (Pn-l,j - Pn-1,/) 717 Wft.-i.lCPn-l.j) 

is right eigenvector of B n (X): 
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n-l 



B n (X) |* P J = A r nfi uj~ pn '° H (A + r n , k u-^ k ) |* Pn > (9) 

fc=i 

if |^ r Pn _ 1 ) is right eigenvector of _B„_i(A). The Fermat curve points p n , p n j, p n j, p^, l and 
r n fc entering ([HD, (JHD are related to the parameters of the model a s ,b s ,c s ,d s , k s by equa- 
tions (54)-(61) of [1], e.g. p n ,i = (x n ,i,yn,i), x n ,i = d n /(x n c n r nil ), = r„ )fc /r n _ M . 
For the more involved determination of the r^j, see Section 2.3 of [1]. 
Analogously, the vector 

(*pJ = J] Q L (Pn-l,Pn|pn)(*p„_ 1 | ® n(^J (10) 

where 

n-2 

Q L (Pn-l,Pn|Pn) =^- ( ^-^- l)(P "-"-° )+ ^- 1 -^- 1 ' TT IU p »-l,l (p n _ij - p n -l,, - 1) X 

3,1=1 

«Vo (Pnfi ~ Pn-l,o)Wp n (Pn ~ Pn) "TT (Pn-1,*) 



X 



f-n , ,^ ' ^ ^ 



n-l 



nr=i nLi (pn-^ - p«,fc - 1) ^ ^-^(-pn-i,*, 

is left eigenvector of B n (X): 

(* Pn \B n (\) = Xr nfl uj- p "' J] (A + r„, fc u;-^) (* Pn | 

fc=i 

if (\I/p n _ 1 | is left eigenvector of _B n _i(A). The definition of Q L (p n _i, p n |Pn) uses the same 
Fermat curve points p n , p n j, p n j, pffi as were used for the right eigenvectors. The proof 
of ( fTTT) by induction is lengthy, but analogous to the proof of (jHJ) given in pQ. 

At A equal to one of the n — 1 zeros of the eigenvalue polynomial of B n (X) 

K,k = -r n ,kU~ pn - k , k = l,...,n-l, (12) 

the operators A n (X njk ) and D n (X njk ) act as shift operators for the fe-th component of 
the vector p n , -D n (A n> &) shifting in addition also the zeroth component: 

An(\n,k)\* Pn ) = MPn) 1^), (*pJ4i (An,*) = UJ^^Pn") (* p "* I , (13) 

D n (\ n , k )\y pn ) = MPn) l*p+°-- fe )' (V Pn \D n (\ nik ) = uy k {p'+ k ) (* p -o,+*| , (14) 
where ~ n ~ 2 

Tn _ s=l 

r uPn-Pnfi-^ 

MP'n) = --fT- n,k [[ F UKk) (15) 

' n - 1 lls=l Vn-1,8 m=l 

and 

F m {X) = (b m +u)a m x m \) (Ac m + d m j ' x m ) . (16) 

The operator V n = V1V2 . . . v n , which defines the term of highest degree in A in A n (X) 
and the free term in D n (X), shifts the zeroth index of vector p n \ 

VnltfpJ = |* p +°>, (*pjV„ = (^ p -o|. (17) 



Form-factors in the BBS model I: Norms and matrix elements 



6 



Using interpolation polynomials and formulas (ITS]) . <HM and (TT7j) . one can construct 
how A n (X) and D n (X) act on the left eigenvectors of B n (X), analogously to what was 
done in equations (66)-(68) of [TJ for the right eigenvectors. 

3.3. The norms of the eigenvectors of B n (X) 
The pairing k (i\l)k = <V, 7 implies (we use y k = = uy%): 

N fx k \ N ' X 

k(ipp>\^ P )k = <V,p — — • (18) 
u WkJ 

Proof. From (j3J) we get directly 

— rr / — tt 

since w^j) = w^i^) tu 7 . States for p' ^ p belong to different eigenvalues in (jHj), (EJ), 
so they are "orthogonal" (have vanishing pairing). We then use ([3]) and 

N N-l 1 , r 

_^ = y ^,7' y 1 = £ 

1 - cjt x ^ ' ^ y N ' 

7'=0 7GZjv * 

and for 1< a < N: V — = N . □ (19) 

Next, we calculate the left-right overlap of the eigenstates of B n (X) for general n. 
We shall consider only the generic case of parameters such that all the eigenvalues of 
B n (X) are different (in particular, this is important for the action of A n (X) and D n (X) 
on the eigenvectors of B n (X)). Therefore (^/ p i|\l/ p fl) = if p^ ^ p^. The following 
theorem gives the value of (^pl\^/ p r) ("norm") for p^ = p^ = p n : 

Theorem 1 The norms (ty pn \ty pn ) are independent of the phase p n fi and their 
dependence on p' n = (p n ,i, ■ ■ ■ , Pn,n-i) is given explicitly as 

C n C Tl 

where the normalizing factor C n is independent of p n and it is defined recursively by 



n (r N - r N ) TT (r N - r N 1 n ~ 2 
n n \T AT n-2 i-i-k<k'\' n,k n,k'J 1 IKl'V n-l,l' n— 1,2/ "1 T N+l / 1 \ 

<~>n ~ O n _i7V n iV TT ( r N _ r N \ ll r n-l,/' 

llfc,A'n,fe ' n-l,l) l=1 

C\ and C*2 are given by 

N ( x i\ N1 „ ^( *2 xN -' 



u \yij uj \y 2 y2V2,Q 

and 

N—l / n— 2 ftii- 1 n,A; 
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Proof. We shall give the proof by establishing an explicit recursion C n _i — > C n . For 
n = 2 there is only a single zero r 2t iU>~ P2 ' 1 and the denominator in (1201) is unity 
Similarly for n = 1 this denominator is unity too as we see from ( fl8l) . 
From (J?D and (fTUj) we get 

(*pj*p„) = ^(Pn-bPnlPn) Q R (Pn-l'P"IPn) (^Pn-ll^Pn-i) n^pJV'pJn- (23) 

P n _ 1 e(z A O n - 1 

Inserting here the explicit expressions (iTTj) for Q L , and jS]) for Q R , using (fl8|) . and 
performing the summations over p n -i,o and p„ by (|T9l) we get: 

n_2 FT"" 2 (V 1 • uj~ Pn - 1 ^ — r MU)~ pn - 1 A 

x n(r „_ M ^-v) 2 ik £ ^ W ,,^,,), (24) 

where we combined several phase-independent factors into the quantity (j22l) . The main 
issue now is to perform the summations over the phases p' n _i explicitly. In order to 
avoid too many indices, let us define 



a 



W=( ai ,...,a n -2) = (r„_i,ia;-^- 1 - l ,...,r ft _ 1 , ft _ 2 w-^- 1 - tt - a ); 

b^=(6i, . . . , 6 n _ x ) = (-A n ,i, . . . , -A n>n _x) = (r ntl u- pn '\...,r n>n ^u- p ^- 1 ). (25) 

Then using the induction assumption that the formula for the norms is given by f[2"Uj) . 
the relation (12^1) reduces to 



C n ^ _ U^-at) CU 



Ui <m (b m - b t ) * A/ " Pn _ l l) 5„_ 1 ,„_ 2 H " ^ 'lTLi( a »- 6 *) n^= 2 /+i(am-aO 

^^KC^^aW.bW), (26) 
where we define f' n _ x = r n _i ; i . . . r n _i jn _ 2 and 

S„(a<»\b<»>) = J2 nU U ^h< a ' I^ ) ■ P7) 

In the Appendix, we prove the following identity 

( n ~ 2 \ n_1 h N h N V\ n ~ 2 (n N -n N ) 

q (An) K (n)\ _ yvrn-2 T~f n N \ TT °fc ~ LLl<l'\ a l a l> ) / 9S n 

The (Z7v) n_2 -summation in (I27p runs over the discrete phases of the variables aS n \ 
From the Z^-symmetry the result of this summation ff28l) can depend only on the iV-th 
powers of these variables. Substituting (T28]) into ff26l) . and using the notation fT25l) . we 
get the recursion ( I2T1) . □ 
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3.4- Action of u„ on B n (X) -eigenvectors 

Now we shall calculate an explicit expression for the action of u n on the B n (X)- 
eigenvectors \^ Pn )- From the relation 

taking into account a n /r n = l/x n , b n /r n = y n /x n , we get the formula for the action of 
u n on one-site eigenvectors: 



U n|VV>n)™ = ^ 



f-|W»--|^ + i>»)- (29) 

We use this formula to obtain the action of u n on general eigenvectors of B n (X), which, 
in particular, will show that most matrix elements of u n between eigenvectors of B n (X) 
vanish: 

Theorem 2 The action of u n on eigenvectors of B n (X) is given by 

un|^) = ^\* Pn ) - u^y*^^) + (30) 



n—1 / ~ ~ \ I ttl — 2 n.k 
+ y^uA.0-& ( VnTn - X nfi r n-XVn \ 1 i/=l Vn-1,1 ^ 



o +k l 

I'll 1 



Proof. To prove ( !30l) we use (JZJ), ( 1291) and rewrite ( l30l) as a relation for Q R (p n ~i, p n \pn)- 
The terms of the resulting relation can be separated into two groups, giving two relations 
which can be proved independently: 

= + E „ , n:L "' 2 ^ q r (p.-.,^i^) (31) 

f-f r„ ifc ^ fti ' fe r„ i0 ll(-r„ )fc a; + r„ iS ^ Q K (p n _i, p n \p n ) 



and 



lQ R (pn-l,Pn - l|pn) p„ „ Q R (pn-l , Pn | p+°) 



U — TvrT i \ _ u Vn.O 



Q R {pn-l, Pn\Pn) ' Q R (Pn-l, Pn\pn) 

+ jjPn-O-P"- 1 Xn '° T "~ 1 11;=1 Vn-1,1 Q R (Pn-l, Pn\Pn k ) ^) 

^ r n _ ll0 EI (- r n,fc^" P "' fc +^n,^~ Pn ' s ) Q R (pn-l,Pn\pn) 



k=l 

sfk 

In order to verify both fl3Tl) and fl32l) we start evaluating the ratios of the Q R using the 
explicit formula (jHJ), e.g. 

since in ( 1331) the shift p+ fc affects only three terms in (jSJ) and we use ([3]) and 
^n-i/ = r n,k/ r n-i,i- For (l3T|) . inserting (1331) . and after some cancellations, collecting 
the fc-dependent terms, the sum over k can be performed using the identity 

E , nH^"^^ = FT^f & = -r^-^, = -r n ^---. (34) 

fe=l lls=l,s^feVSfe SsJ llfc=l Sfc 
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For verifying (I32p we have to calculate the other shifts of Q R too. This time, with the 
same and d as before, we perform the sum over k using the identity 

\p nr=i 2 (6c - CO _ 1 n 

4. Periodic model. Baxter equation 

To treat the periodic model it is convenient to make a Fourier transform over p n fl of the 
eigenvectors of B n {\). This yields a basis of eigenvectors of the operator V n of ( fT7l) : 

(* AA |V„ = w' V„|* ftA > = w>|* ft *>- (35) 

V„ is the total Zjv-spin rotation operator and p G I^n is the corresponding total charge. 
Let ($ p ,e| and |$ p ,e) be left and right eigenvectors of t n (A) with eigenvalue 

t n {\\ p, E) = Eq + E\\ + h ^-iA"- 1 + £ n A n , (36) 

where E = {E±, . . . , E n -±} and the values of E and i? n are 



n 



E = l+^H b -^, E n = f[a m c m +uj?l[x m . (37) 

m=l m m=l m=l 

Following the general procedure of the SoV method, we express these eigenvectors as 
($p,e| = Q L (Pn\p'V) (*p,pd, |$ p ,e) = £ Q R (p' n \p,E)\y p , p ,J, (38) 



wherdjj 



n— 1 n— 1 



g L ( P ;i P ,E) = J] g^( Pri|fc ) J] «, y^-^-i), (39) 



fc=l s,s'=l 



Q R M P, E) = " fc=1 " fcTO/ 7 • (40) 



ilLl Qf(Pn,k) 
n"T'=l W n n,s'(Pn,s - Pn,s<) 



Pn , S 
s^s' ' 



The products of Zjv-cyclic functions w n ,s>(p n .s — Pns') hi (13"9"j) and (HO]) are Sklyanin's 

Pn,s 

measure which makes the rest of the kernels Q h (p' n \ p, E) and Q R (p' n \ p, E) factorizable 
into products of single variable functions (SoV). The Baxter equations for these functions 
of separated variables are 

t n (X n>k \p, E) Q\{p n ,k) = u n ~ 1 A k (A n , fe ) Q\{p n ,k + 1) + c^A+^A^) Ql(p n ,k - 1) , (41) 
t n (X n , k \p, E) Q^(p„, fe ) = A+(A n , fe ) Q R (p n , k + 1) + A* (wA*,*) Q^(p„, fe - 1) , (42) 

| The kernels in {35]), (]5§j) and (HP)) are not the kernels (O and (jlip used in the recursive definitions 
of the left and right eigenvectors of B n (X), observe the different types of arguments. 
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where X n ^ = —r n ^w p7l,k , k = 1, . . . , n — 1, are the zeros of the eigenvalue polynomial 
of B n (X), see (JT2D, and 

n-l 

A+(A) = (w'/xk) (A/o;) 1 -" J] F m (A/w) , A fc (A) = Xfe (A/^ 1 F B (A/w) , 

m=l 

n-l n— 2 

r n,0 r n-l /TT ™,s / n.k\ TT "ife /in\ 

Xfc = — — (11 ynlk/Vnis) 11 y n -i, s - (43) 

8=1 8 =1 

5. Periodic BBS model for N = 2 

5.1. Eigenvalues of transfer matrix for generalized homogeneous Ising model 

In this section we consider in more detail the case of the N = 2 periodic homogeneous 
BBS model, where u = — 1. By homogenous we mean that the parameters a, b, c, d 
and h each are taken to be independent of the site index. As it was shown in [10J this 
model is a particular case ("free fermion point") of the generalized Ising model. 

In the N = 2 case, we have u k = a% and v fc = erf, where cr|| and a% are Pauli 
matrices acting on the spin at the fc-th site. The operator V n = a\a\ ■ ■ -a^ is the 
spin-flip operator. It commutes with the transfer-matrix t n (A), and we have simply 

= 1. Therefore the eigenvectors of t n (A) are divided into two sectors according to 
the eigenvalue (— l) p , p G {0, 1} of V n . If the parameters a, b, c, d and x are generic, 
no degeneracies of the eigenvalues occur. The transfer-matrix of the standard Ising 
model on the finite lattice can be obtained from the N = 2 BBS integrals of motion 
Hfc, k — 0, 1, . . . , n, when parameters of the homogeneous N = 2 BBS model satisfy the 
relations 

a = c, d = b, x=l and X = b/a (44) 

(see Section 16. ip . In this case degenerations of the spectra occur and in order to 
distinguish the corresponding eigenvectors one has to use an additional operator which 
commutes with transfer matrix. In case of the periodic model, such an operator is the 
operator of translation by one site: 

T |7l)l ® |7 2 ) 2 ® • " " ® \jn)n = |7n)l ® |7l)2 ® • • " ® |7n-l)n , 

T ^ T- 1 = [ T, t n (A)] =0, [ T, V n ] = . 

Our construction of the eigenvectors fl38|) is obviously non-invariant with respect to 
translations and we were not able to show that the vectors ( 1381) are also eigenvectors 
of T. A translational invariant description of the eigenvectors for the generalized Ising 
model has been given in [18] . The description of the spectra of the transfer matrix given 
in this paper coincides with the description of the eigenvalues obtained in our formalism 
presented in [TJ. This allows to identify the eigenvectors in both descriptions. In 
particular, the eigenvectors are labelled by quasi- momenta of the excitations q = ±^m , 
where m is integer or half-integer depending the eigenvalue (— l) p of V n . So the whole 
space of states decomposes into two sectors according to the value of p: 



Form-factors in the BBS model I: Norms and matrix elements 



11 



• NS-sector: p = 0, the eigenstates of t n have an even number of excitations with 
quasi-momenta q G NS = {^(Z„ + 1/2)}. The dimension of this sector is 2 n_1 . 

• R-sector: p — 1, the eigenstates of t n have an odd number of excitations with 
quasi-momenta q G R = {— Z^}. The dimension of this sector is also 2 n_1 . 

In [T] the eigenvalues t(X) (we suppress the chain-length index n) were derived from a 
functional equation which could be written in the form (don't confuse these A(q) etc. 
with the A n (X) etc. in (ED): 

t(\)t(-x) = (-irn(^(q)^ 2 -c(q) +2»*As( q )), (45) 

q 

where ^( q ) = fl 2 c 2 _ 2x ac cos q + x 2 , fi(q) = (ad - be) sinq , 

(46) 

C(q) = 1-2(6 d/x) cos q + b 2 d 2 /x 2 . 

Solving (1451) . the 2 n eigenvalues corresponding to the NS- and R-sectors are given by 
(J4"7j) [U [TOl [18] where we have to choose all possible sets of ±— signs (which we write 
(-l) 47 " with a q G {0,1}): 

NS: t(X) = (a"c" + x«)n q6NS (A+ (-l) CT ^ q ), n^NS^ 1 )^ = + 1 ' 
R: t(X) = (aV-x«)n q6R (A + (-lN, IlqgR (-1)" = ~1, 

the restriction on the signs follows from fl35l) . We shall also write A q = (— l) CTq s q , and 

the amplitudes s q are given by 

s q = (VS(q) -iB(q))/A(q), D(q) = A(q) C(q) - 5(q) 2 . (48) 

The convention which fixes the sign of ^ D(q) is given in (120), (121) of pQ. States 
which have eigenvalues where in ( 1471) there are some minus signs (— l) CTq = —1 are said 
to have excitations with quasi-momenta q. The restriction on the signs means that for 
the NS (R) sector the number of excitations, i.e. the number of minus signs, must be 
even (odd). 

As an example, the 2 3 eigenstates and eigenvalues for the 3-site chain (n = 3) are: 

• NS-sector (p = 0): 

|)ns ( no excitations) — > t(X) = (a 3 c 3 + x 3 )(A + s 7r / 3 )(A + s^ 7T / 3 )(X + s w ), 
| ± tt/3, 7r) NS , -> t(X) = (a 3 c 3 + x 3 )(A =F s 7r/3 )(A ± s_ 7r/3 )(A - s n ), 

|tt/3, -7t/3) ns , -> t(A) = (a 3 c 3 + x 3 )(A-s 7r/3 )(A-s W3 )(A + s 7r ). 

• R-sector (p = 1): 

|0) R -> t(A) = (a 3 c 3 -x 3 )(A + s 27r/3 )(A + s_ 2 ^ /3 )(A-s ), 

|±2tt/3) r -> t(A) = (a 3 c 3 -x 3 )(A T ^ /3 )(A± S _ 2V3 )(A + So ), 

|0,2tt/3, -2tt/3) r -> t(A) = (a 3 c 3 - x 3 )(A - ^ /3 )(A - S _ 2 , /3 )(A - s ). 

From [18J, it follows that the action of operator of translation T on the eigenstates 
gives eigenvalues of the form e lP , where P is the sum of the quasi-momenta of all 
excitations of the given state. In case of the Ising model when ad = be and B(q) =0 we 
have s q = s_ q : two states which have excitations with quasi-momenta q and — q have 
the same eigenvalues but can be distinguished by the eigenvalue of T. 
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The eigenvalues of t n (A) provide the existence of nontrivial solutions of the systems 
f HTj) . f|42|) of homogeneous Baxter equations. These solutions give explicit formulas fl38l) 
for the eigenvectors, which form a basis in the space of states of the periodic BBS model 
for N = 2. 



5.2. Norms and orthogonality of the eigenvectors of the periodic N = 2 BBS model 

Let us fix an eigenvalue polynomial t(A) of A n (X) + D n (X) corresponding to a right 
eigenvector from the sector p. In order to find this eigenvector explicitly we have to 
solve the associated n — 1 systems (k = 1, 2, . . . , n — 1) of Baxter equations: 

t(-r n , k )Q*(0) = (A+(-r n , fe ) + A-(r n , fc ))Q*(l), 

(49) 

«W «(1) = (A+(r n , fc ) + A-(-r n , fc ))Q^(0). 

Since t(X) is eigenvalue polynomial, it satisfies the functional relation. In pQ it was 
shown, that the functional relation ()45l) ensures the existence of non-trivial solutions to 
( 1491) with respect to the unknown variables (0) and Qf(l) for every k — 1, 2, ... , n—1. 
In the N = 2 case, this means that for every k we have one (in the case of 
degenerate eigenvalues, possibly zero) independent linear equation. In the case of generic 
parameters, both hand sides of each equation are non-zero. So we may fix Qf(0) 
and obtain two equivalent expressions for Qf(l): 

qR,^ = t(- r n,k) = Afc(w) + Afe(-r n , fc ) 

A+(-r n , fe ) + A^(r„ ifc ) t(r n>k ) 

Similarly, for the left eigenvector we have 

t(-r n , k )Ql(0) = (-ir- 1 (A+(r„, fc ) + A,T(-r n , fe ))QL( 1); 
t(r n , k )Ql(l) = (-ir- 1 (A+(-r n!fc ) + A,T(r n , fe ))QL( ). 
Fixing Q k (0) = 1 we obtain 

G L m (-l)»-it(-r n , fc ) A+(-r w , fc ) + A fc -(r n , fc ) 

) A+(r n , fc )+A fc (-r nife ) (-1)^-1^) • 

Since for generic parameter values we shall have t(r nyk ) ^ 0, these explicit formulas give 
immediately 

Q L k (Pn, k )Q^(Pn, k ) = (-lr^-^u-ir^r^AK,) . 

Eigenvectors in the periodic case are defined by the formula In the case N = 2 
we have a simple inversion relation for the cyclic function w p (p): w p (p)w p (p — 1) = 
y/ (1 + x). This allows us to use instead of (|59l the following kernel 

Q\Pn\ P, E) = QliPn ' k) , • (51) 



(50) 



EI™ s'=l W n V:' s ' (Pn,s ~ Pn,s') 



§ When the parameters satisfy the Ising model restrictions, it is not always possible to choose such 
normalization. 
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Up to a coefficient, (IBTj) with (I3"5j) gives the same eigenvectors of A n (\) + D n (\) as 
(139|) . In what follows we shall use only fl5Tl) . The pairing of ($| and |$), which are 
left and right eigenvectors corresponding to the eigenvalue t(X), gives (we will use the 
normalization by (^0,0 1 ^0,0) to get rid of unimportant factors appearing in the formula 
for the norms) 



E 



<*0,o|*0,o) ^ n?<m( W P??(PnJ ~ Pn,m)w p % l m (Pn,m ~ Pn,l)) 2 <*0,o|^0,o) 

where for the norm we have used (1201) which leads to 



1,1 



rn-l/ \ ' 



where = (0, 0, . . . , 0) has n — 1 components. In the N = 2 case we have 

y 1 — x 1 1 + x 



(52) 



w p (l) 



1 + x y ' (w p (l)) 2 1 — x 



Therefore for the Fermat point p™'™ = (a^'™' 2/™'™) with coordinate x^™ — r n,m/ r n,i we 



Wv n ' m {Pnl- Pnm)) = (~1) ' " ! 7 — TT 



and so 



2 



-1) 



Combining all these formulas we get for the left-right overlap of the transfer matrix 
eigenvectors of the periodic BBS model at N = 2: 

(*o,oi*o,o> n£*(»w) ^ nr<™((-i) p -^ + (-i)^^)' 

Using the same techniques as we used for calculating the norm in the auxiliary 
problem, we are now able to perform the summations in ( 1531) explicitly: We write the 
polynomial t(X), equations (I36I) . (14"T|) . as 

n 

*(A) = A[](Hft), (54) 
fc=i 

where A = a n c n + (-l) p x n and p k = A q with q = ^{k + 1 ^ E ). Then the sum in ([53]) is 
A n-1 JJr y n , llLl((-) Pn -' r ^ + Pk) 



i=i i 

11 rr™- 1 ^2 _ ^2 \ -Z^ll' 



TT llA:=l( r n,f /4) 2 „_l -/ 1 T f4 ~ ^ lL</' ( r nJ r n,l') 



A n-i 2 „-i ^ "Q ^ + ^ ? = rn XTn 2 . . . rn n _ 1 . (55) 
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In the second line the sum takes just the form (j79p which we had for S n+ i with N = 2 
and a( n+1 ) = (r^-)**- 1 , . . . , r fl>T ,_i(-) ft *-»- 1 ) and t>( n+1 ) = . . . ,//„). In the last two 
lines we have used the result (IHUj) for S^+i derived in the Appendix. So we obtain the 
norm of the general periodic state vector 

T^TT = 2- 1 * n F < ^l-r + r :' ) , f[C* + H) ■ (56) 

<*0,o|*0,o) \Ak=l\Al=l ( r n,l + M X<j 

If the spectrum of the model is not degenerate, the overlap of the eigenvectors 
corresponding to different eigenvalues is zero. In case of degeneration, as it happens in 
the case of the parameterization a = c, b = d and x = 1, corresponding to the Ising 
model, one should be careful using (1561) . In Section [6] this case will be discussed in detail. 

5. 3. Matrix elements between eigenvectors of the periodic N = 2 BBS model 

In Theorem 2 we found the action of u n on an eigenvector |^/ p „) of B n (X): the result is 
a linear combination of the original vector plus a sum of vectors which each have one 
component of p n raised. In order to get the matrix elements of u n in the periodic model, 
we first obtain the matrix elements between Fourier transformed eigenvectors of B n (X) 
defined in f l35|) 

<*ai<J = <*o X | + (-r<*ixl , I***) = |*o,A> + H P l*iX> • 
Using ( l30i) . the following n-site matrix elements of u n are non-zero 

^ n I 1 _L ^ n ' C ^ n _ a n , -t\p' , x l x 2 ' ' ' K n-\b n 





u„| 




(*0,p{J 


Un| 


r o,p;) 



_r_XjPn _|_ — J!/ nPn _|_ 

•"n 1 n 1 n,0 



Xn ( — l)Pn _ Un^y™ _ an (_-[\P'n _ Xl><2 ' ' ' X n-lb n 



JVn.O V x nC n r n ,fc / n^fc( r n,fc( _1 ) Pn: ' i + ^n,s(-l) Pn ' s ) ' 

where p' n = p n> i + p riy 2 + ■ — h p n ,n-i- Note, since u n anti-commutes with V„, all matrix 
elements of u n between the vectors from the same sector p are zero. 

Then by (138]) we pass to periodic eigenstates transforming by the solutions Q of 
the Baxter equations. Let (<&o| be a left eigenvector of the transfer-matrix t n (A) with 
p = and |$i) be a right eigenvector with p — 1. Let (/0 n ,fc) and Q% (p n ,k) be 
the solutions of Baxter equation corresponding to these two eigenvectors. Then for the 
matrix element after some simplification we have 

($ |u n |$i) = V^/_ 1 x n ^ TT r n,l + r n ,m y 

(*o,o|*o,o> ~ 4" iJn (- 1 ) Pn ' lr n,i + (-^ m r n , m 
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x 



7LT \ Tn r ™>° / 



1=1 

n— 1 n— 1 



jt=i i^fc 

x an&nC w A | (-l)'-*rf w \ 1 \ (5g) 

r «,0 V X nC n r n ,fc / ris^fc( r n,fc(- 1 ) Pn '' : - r n,s(-l) Pn ' s )/ ' 

The product on the right hand side of the first line comes from the change of 
normalization from (jSTj) to fl58|) by the factor (^o.p'l^o.p')/ (^o,o|^o,o), which can simply 
be read off from ( l20i) with C n cancelling. 

We have not yet been able to perform the summation over p' n in (|58|) for general 
parameters of the BBS-model. However, for the homogenous Ising case (1441) . we can 
show that (1581) can be put into a fully factorized form, although in this case special 
complications appear from the coincidence of the zeros of the transfer matrix with the 
polynomial zeros of B n (X). In the next section we give the result of this summation, 
which proves a formula for the matrix elements of the spin operator (form-factors) on a 
finite lattice conjectured [HI [15] by A. Bugrij and O. Lisovyy. Details of the proof which 
is a new result for the Ising model, and the comparison with the notation of [HI EE] are 
relegated to the sequel article [16J. 



6. The homogeneous Ising model 

6.1. Relation to the standard Ising model 

In this section we restrict the parameters of the homogeneous N = 2 BSS model to be 
a = c, b = d and x = 1. So the cyclic L-operator (CQ) reduces to 

L U {\)=( \ +X l\ X Ti a ~ b l k) )- (59) 
V ; ^ u fc (a -6v fc ) \a 2 + v k b 2 J v ; 

Let us make the special choice of the spectral parameter A = b/a as in (jHJ). Then the 
L-operator degenerates and we get 

L h (b/a) = (l+v k b/a)( 1 bU " ) = (1 +v k b/a) ( 1 )(l bu k ) . (60) 
At this point the transfer matrix is 

n n 

t n {b/a) = tr L x (bla)L 2 (bja) ■ ■ ■ L n {b/a) = JJ(1 + v fc • b/a) ■ JJ(1 + u fc u fc+1 -ah). (61) 

k=l k=l 

Recall that due to the periodic boundary conditions u n+k = u^. Using 
exp(i^ x u fc u fc+1 ) = coshiC a .(l+ UfeUfe+itanhifj.), exp(i^*v fc ) = coshK*(l+ v fc tanhK*) 
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and writing u k = a z k and v fc = crjF, it is easy to identify t n (b/a) with the standard 
Ising transfer-matrix: 



tising = exp ( K a k ) ex P ( Yl Kx a k a 

b 



k+l 

k=l / \k=l 



e 



- 2R y = tanh K * = - , tanh K T = ab . 



a 



where and K* y are coupling constants of Ising and dual Ising models on the square 
lattice along the X and Y axes, respectively. 

The eigenvectors are obtained by the method of separation of variables, they do 
not depend on A. In what follows we shall not fix A = b/a and so we consider a family 
of models depending on A, a and b. This family includes the Ising model at A = b/a 
and the eigenvectors obtained by SoV are eigenvectors of the Ising transfer matrix too. 
Let us note that this L-operator formulation can be extended [10] to a larger family of 
Ising-like models, giving a possibility to use SoV for finding explicit formulas for the 
eigenvectors of corresponding transfer-matrices. 

In ( )45l) and ( )47l) from [lj we already quoted the functional equation and the 
eigenvalues of the general homogeneous N = 2 BBS model. Introducing 

T(A 2 ) = 1 + b 2 d 2 /x 2 - A 2 (x 2 + a 2 c 2 ), F(X) = (b - axA) (Ac + d/x) , (62) 

the equation A(q)X 2 — C(q) + 2iX q B(q) = can be written as 

T(s 2 q ) = e^F(s q ) +e ^F(- Sq ). (63) 

The solution of (1631) is given by (H8l) where q = ±- m, with m even for the R-sector 
(p = 1) and m odd for the NS-sector (p = 0). 

The equations for the amplitudes r n ^ of the roots of the eigenvalue polynomial of 
the operator B n (X) of equation dHJ) were derived in equation (A7) of [ljjjj] 

T ( r l,k) 2 = 4cos 2 (g nifc )F(r nifc )F(-r nifc ) , g n>fe = vrfc/n, fe = 1, 2, . . . , n - 1. (64) 

In the case of the Ising model parametrization (144p we have -F(A) = b 2 — a 2 A 2 , and 
([63]) and (JMD reduce to T(s 2 ) = 2cos(q)F(s q ) and T{r 2 nk ) 2 = 4 cos 2 (g n , fe ) F 2 (r n , fc ) , 
respectively. Then the solutions of these equations are: 



fr 4 -2fr 2 cosq + 1 b*-2b 2 cosg n , fc + 1 

Sq ^a 4 -2a 2 cosq + l' ^ y a 4 - 2a 2 cosg n , fc + 1 ' 1 j 

where the momentum q in both NS- and R-sectors takes the values {0, it, ±g n ,&} with 
< q n ,k < 7r, in particular g n> fc 7^ 0, n. So we have 

6 2 -l b 2 + l 

r n ,k = s qn k = s. qn k , s = - , s w = ^q-y . (66) 

The possible coincidence in the Ising case of the two quite different parameters: the 
zeros r n .k of the polynomial B n (X) and the zeros s q of transfer matrix will create some 
peculiarities when in the following we set out to solve the Baxter equations. 



In p] we used n ,fe/2 instead of q n 
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6.2. Solution of the Baxter equations 

We have seen that in the Ising case ( )44l) the eigenvectors and eigenvalues of t n (A) 
decompose into two sectors p = 0, 1. With (jSJ) in addition we have also F(X) = F(—X) 
and the Baxter equations (HT|) . (H2|) for left and right kernels become identical. Omitting 
the superscripts L and R on and writing X U)k = — (—l) Pn ' k r ntk , p n ,k = 0, 1 we obtain: 

KK,k) Qk(Pn. k ) = ( { ~!l PF l + (-An^)"- 1 Xk F(X n , k ) ) Q k ( Pn>k + 1). (67) 

In order to solve these Baxter equations for a given sector p, we need t(±r nyk ). Now 
from (H71) we have 

t(A) = (a 2 " + (-1)") n q (A + A q ) , A q = (-ip Sq , (68) 

and we see that i n (±r n> &) vanishes if A q = ~^r nk . When p and & have the same 
parity (— ) p = (—) h the quasi-momenta q which describe the spectrum A q of the transfer 
matrix do not coincide with any q n ^ k and the value of the transfer matrix t(±r n>k ) does 
not vanish. If p and k have different parities: (— ) p = (— ) k+1 this could happen. Below, 
when we solve the Baxter equations we shall treat these cases separately. 

Moreover, as we can see from fl66l) . for the Ising parameterization the amplitudes 
s qnk = S- Qnk coincide. This leads to a degeneracy of eigenvalues of t n (A): an eigenstate 
with q = q n> k excited and q = — q n ^ not excited has the same eigenvalues as the 
eigenstate with — q n ^ excited and q n ^ not excited: both have the same factor (A 2 — Sg n fc ) 
in the eigenvalue polynomial. The degeneracy of eigenvalues of the transfer-matrix 
can lead potentially to the following problem: the functional relation guarantees, of 
course, the existence of a non-trivial solution of the Baxter equations for any particular 
eigenvalue-polynomial, but we need in our case two independent solutions. Fortunately, 
due to the coincidence s q= ± qn k = r n fc both sides of the Baxter equations for the unknowns 
Qfc(0) and QkiX) become zero and we can build a two-dimensional solutions space. 

The compatibility condition following from fIBTl) is 



t(-r n , k )t(r nik ) = (-1)"- 1 ( V + (-r n , k ) n - l XkF(r n , k 

where t(A) is an eigenvalue from the sector p. If k is such that (— l) h = (— l) p+1 , 
then the quasi-momentum q = q n ^ k belongs to the sector p and for r n ^ = s qn k we have 
t(— r ni fc)t(r n> fc) = 0. This implies a relation not depending on a particular t(A) and its p: 

xlrZ- l) = ^Y +k+1 F n -\r nM ). (69) 

In order to find the eigenvector corresponding to t(X) from the sector p we have to 
find all the Q k (p n ,k) for all k solving the Baxter equations. We need to distinguish the 
following four cases with respect to the value k: 

(i) (— l) p = (— l) fe . In this case both t{r n)k ) ^ and t(— r n>k ) ^ 0. So we may fix 
Qfe(0) = 1 and, using ( 1691) we obtain 

(_l)«-i t (_ r ) 



n— 1^ • '•■ - 

n— 1 



Qk(l) 



SXfcC 1 F(r n , k ) 
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The other three cases correspond to (— l) p = (— l) k 1 so that the big brackets on the 
right-hand sides of the Baxter equations are zero due to (1691) . 

(ii) t(r nj k) 7^ 0, t(—r njk ) = 0: The eigenvalue polynomial £(A) contains the factor 
(A + r n fc) 2 (both momenta q = iq n k are not excited, i.e. both not in the spectrum) 
and we have 

Q fc (0) = 1, Q fc (l) = 0. 

(hi) t{r n)k ) = 0, t(— r n>fc ) ^ 0: t(X) contains the factor (A — r„ fc ) 2 (both momenta 
q = ±g„ fc are excited) and we havejjl 

Qfc(0) = 0, Q k (l) = 1. 

(iv) both t(r n> fc) = and t(—r n ^) = 0. This case happens when only one of the two 
quasi-momenta q = ±q n , k is excited, so that t(A) contains a factor (A 2 — r 2 k ). As we 
explained above, this eigenvalue is degenerate. In this case both hand sides of both 
Baxter equations are zero. In principle, we can choose any two independent solutions 
for Baxter equations, but, in general, they will not give us eigenvectors of operator of 
translation T. Since our construction of the eigenvectors ( 1551) is obviously non-invariant 
with respect to translations it is unclear how using the direct action of T one can choose 
solutions of the Baxter equations which give such eigenvectors. We take another way: 
In order to obtain eigenvectors of transfer matrix which are eigenvectors also of the 
translation operator, we shall first lift the degeneracy by starting with parameters such 
that ad — be = rj (keeping x = 1) with t] small but finite. At the end we shall take the 
limit 7] —>■ 0. Observe that for rj ^ (|49p and (150|) are different. From (j49H we get 

O r (()) = 1 o R m = t(- r n,k) 

^ k 1 ' ' Wk 1 ) (-l)VOa Kfc)"- 1 ) F^(r n , k ) + (-l)n-i Xfc ( rn>fe )n-i F(-r n , k ) ' 

When rj — > 0, both numerator and denominator of Qj^(l) tend to 0. Let us find their 
leading (in fact, linear) terms in rj — > 0. For the denominator we have for 77 — > 

' —l) p F (r n ,k) , j , / \n~l 771/ _ \ „ / i\n+l^_ ,. /_ \n-l 



1 



Xk(r n ,k) n 1 

where we used the relation (valid for arbitrary rj) 

F n -\r n>k )F-\-r n>k )= X irZ- l) . 
Let us write the transfer-matrix eigenvalue polynomial ( 1651) as 

*(A) = ^„(A) (A-s qh )(A + s_ q J, 
where q^ = (— l) <T ' !n ' fc+1 <7n,fc is the value of the excited quasi-momentum, so that 

t{-r n , k ) = tq k (-r n ,k) {r n ,k - S-$ k )(r n , k + s q J . 
Using A(q)(X - s_ q )(A + s q ) = A(q) A 2 - C(q) - 2i73(q)A, B(q) = 77 sinq and 
4(3*) r «,fc - C ^k) ~ 0(77) (due to gSD and gHD at 77 -> 0) we get 

t{-r n>k ) rj qfe (-r„ >fc )2i J B(q fc ) r n>fe /A(q fc ) = -t qfe (-r nifc )2isin(q A .) r njfc ?7/A(q fc ) 

In this case we can't normalize Qfe(0) = 1. 
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at 77 — > , where A(q) = a 4 — 2a 2 cosq + 1. Finally, 

qR( ) = 1 gRQ) = (-l) n+f7g - fc+1 2i sin(g n|fc ) ^ (-r n , fc ) 

Using a similar limiting procedure for the Baxter equation (loTTj) we get 

gi(o) = i, gi(i) = -gj(i) . 

Although the Ising Baxter equations were the same for Q L and Q R , in this case, since 

>£(l)andQ£ 



we take the limit of the general non-hermitian BBS case, Qi (1) and (1) come out 



different. 

6.3. Norms and factorized matrix elements for the homogeneous Ising model 

In this Subsection we state a new factorized formula for the matrix element of the spin 
operator in the Ising model on a finite lattice. Let |$o) an d be two eigenvectors 
of the periodic Ising model from the sectors p = and p = 1, respectively. First let us 
consider the case when these states do not contain fermionic excitations with momenta 
q = Qn,k and q = —q n ,k simultaneously, so that the eigenvalue polynomials do not contain 
a factor (A — ry^) 2 . We introduce the subset V C {1, 2, . . . , n — 1} of indices k for which 
the eigenvalue polynomials of |$o) an d l^i) contain the factor (A 2 — r 2 k ), i.e. for which 
we have the case (iv) of Subsection 16.21 Denote by \V\ the size of the set V. 

Let 8 = 1 if o"o = o~ n and 5 = otherwise. We also use the short notations: 
A = (— l) CTO s , X n = (— l) cr ' r '5 7r and = {—iy qn ' k+k q n ,k- Now the matrix element is 
given by the factorized formula: 

($oK|$i) ($i|u n |$ ) 



(^o,o|^o,o) 2 



(A 2 - A 2 pl- 5 )/ 2 (A + K) 5 x 

x tt 2r nfc -r-r r n) fc+r njTO -r-r sin|(gfc — q m ) 

4^( r n,fc + A o)(^n,fc + A^) ^ A r nM -r n ^ m A A sin|(g fc + g m )' 

k,m£T> k,m£T> 

The proof of ( 1701) is given in our sequel article [16] . 

In order to compare (1701) to the results obtained by A. Bugrij and O. Lisovyy [T^fTo] 
we have to calculate instead of ( |7D1 the ratio 

(71) 

which has the advantage to be independent of the particular normalization of |$o) an d 
|<3>i). To do this we have to divide ( JT01 by 

($ |$ )($ 1 |<l> 1 )/(^ 0i o|^o,o) 2 . (72) 

The expression for (1721 was obtained in (1561) for generic parameters. In the case of 
the Ising model the norms of the degenerated states cannot be directly obtained from 
( |56i) . because the numerator and denominator contain zeros. When degenerated states 
appear, that is when T> is not empty, we have to be careful in using the formula (|56|) for 
the norms and use l'Hopital's rule. In order to implement l'Hopital's rule we proceed 



($0! 








$0) 


($o| 


$ )($1 
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similarly to what we did in case (iv) of Section 6.2 for finding the solutions of Baxter 
equations: we go off the Ising point taking ad — be = rj to be small finite and keep 
terms linear in rj. Let qi = (— l) a9n - l+1 q n ,i for / G T>. Then the coefficient of A in 

A(q,)(A -s q! )(A + s_ q J = A(q,)A 2 - C(q z ) +2i£(q,)A (73) 

gives -s q! + s_ q! = 2iB(q l )/A(q l ). Due to A(q,) r 2 n l - C(qj) ~ 0(77) (when rj -»• 0) the 
formula ([73]) at A = r„ jZ gives (r n>i - s q J(A + s_ qi ) ~ 2i.B(qz)r Tli z/ / A(qj). Therefore 

7 — r - — for 77 - 

for the corresponding factor in the formula (l56|) for norm. Note that this result is 
independent of which of the two quasi-momenta q ni i or —q n ,i is excited. Now we have 
to take the product of one term fl56|) for R with another for NS. Since the terms at s^ 
appear only in NS (R) for n odd (even), we get two slightly different formulas for these 
cases: 

For n odd we have 

Ol^l TTVo ^ fife-odd^ r n,k) rifc-even(^0 ± r n,k) 
Z II \ Zr n,k) ■ T=f 77 j C • 7\ 1 7 

Uk<i,k,i'Odd(( r n,k + r n j){±r n , k ± r n ,i)) Uk<i,k,i-e V en(( r n,k + r n ,i)(±r n , k ± r n M 
x * i * , L , (74) 

rife-odd,/-even( V^ :r n,k + r n,l){ r n,k ± T n j) J 



($o| 


$0) ($ll 


$1) 


(*o,ol 


^0,0 


) 2 



and for n even 

n-1 



($o| 


$0) 


$1) 


(*o,ol 


^0,0 


to 



2^n(2r„, fc )-(A + A.)^ 
fc=i ^ 



nfc-cvcn( A o±r„, fc )(A^±r niA 



fe— odd 



(A + r n , fe )(A 7r + 



X 



Ilk<l,k,l~odd[( r n,k + r n j)(±r njk ± r n> /) J rifc<Z,fc,Z-even( ( r «,fc + r n,i) (=t r n,A: ± ?"Vg 

x * 1 * , L , (75) 

rife-odd,/-even( \^ zr n,k + r n,l){ r n,k ± r n> /) J 

where the sign at ±r n>m is fixed '— ' if m e D and '+' otherwise. 

The final result (1711) will be given in terms of A , A^ and r n>k , fc = 1, 2, . . . , n — 1. 
Now if the eigenvalue polynomials contain the factors (A — r n ^) 2 instead of (A + r n ^) 2 
for some I we have to replace r n j — > — r n> j in the final formula for the matrix element 
flUD for all such I. 

Let us finish this section by giving as an example the matrix element for the 3-site 
chain between |<£> ) = |vt,7t/3)ns and |$i) = |0)r. Using (1741) we have 

8r 3 ir 3 2(r 3 i + s n )(r 32 - s ) 



($o| 


$0) 


$1) 


(*o,o| 


^0,0 


to 



( r h - r li)( r 3i - s )(r 32 
Since in this case X> = {1}, formula (170]) gives 

($ | u 3 |$i) ($i| u 3 |$ ) 2r 3 i(s + s n ) 



(^o,o|^o,o) 2 ( r 3i -s )(r 3 i-O ' 
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Finally, 



($ |u 3 |$i)($ 1 |u 3 |$ ) _ (SO + S 7r)( S 27r/3 - sl /3 )(s 2n /3 ~ S n ) 
($ |$0)($ 1 |$ 1 ) 4s 27r/3 (s 2 _ S 2)(s 2 ^/ 3 - S ) 



where due to (16"5|) 




31 ~~ 



6 4 - b 2 + 1 
a 4 - a 2 + 1 ' 



2tt/3 



32 — 



b 4 + b 2 + 1 

a 4 + a 2 + 1 ' 



7. Conclusions 

In this paper we continue our calculation of state vectors and matrix elements of the 
finite-size inhomogenous Zjv-Baxter-Bazhanov-Stroganov lattice spin model, using the 
method of Separation of Variables. In a previous paper [1] we gave the right eigenvectors 
for the auxiliary system and the Baxter equations which determine the periodic 
boundary condition BBS-eigenvectors. Here we complete this work by calculating 
explicit formulae for the corresponding left eigenvectors and norms and matrix elements 
of a single operator. A main result is the expression of Theorem 1 for the norm of the 
state vectors of the auxiliary system: the multiple summation over the intermediate 
states is performed, so that the norm is put into a factorized form. For N = 2 the 
Baxter equations are solved explicitly, and also the norm of the periodic model is put 
into a compact factorized form. Since it is an open task to perform the intermediate 
spin summations for matrix elements of the inhomogenous periodic model, in the last 
Section we proceeded with the homogenous Ising model parameters. There we manage 
to perform the summations and we present a new factorized formula for the matrix 
elements of spin operator, delegating the somewhat lengthy derivation and comparison 
with the corresponding formula in [14j [15] to our sequel paper [16] . 
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Appendix: The main summation formula 

Consider a degree n polynomial f(t) = f n t n + . . . + f and its interpolation formula 
through n + l arbitrary not coinciding points Ci, . . . , c n+ %: 




(76) 
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We shall use ( T75|) for f n = 0, i.e. for polynomials /(£) of degree less than n, obtaining 
a sum rule valid for an arbitrary choice of the parameters {ci, . . . , c n+ i}: 

y f(ck) = (7?) 

Now take an arbitrary polynomial f(t) of degree less than n + N and choose n+N points 
(ci, . . . , c n+ Ar) = (61, . . . , b n , r, ur, . . . , u; 7V ~ 1 r) . Then the sum rule corresponding to 
(EZD is 

since Y[iLo 1 i^ P ( u}Pr ~ u ° l ' r ) = Nu~ p r^ -1 . 

The goal of this Appendix is to prove equation ( 1251) of the main text, i.e. we shall 
show by induction that for the sets ay 1 ' = {ai, . . . , a n _ 2 } = {ria;~ / ' 1 , . . . , r n _ 2 ^~ Pn ~ 2 } 
and b^ n ) = {61, . . . , b n -i} the Z^-symmetrical sum over all phases of the parameters ai in 



*.(.w= e nf ^'^'n ) (79) 

can be performed explicitly, so that we obtain (179|) in the factorized form 

This assertion is correct for n = 3, as can be seen by explicit summation, using 
E P&N b/(ru-»-b) = Nb N /(r N -b N ): 

S 3 (ai] h, b 2 ) = ( ai _ 6i )( ai _ 6a ) = J2 \^h 2 (h-r^-Pi ~ 0l - nu -Pi 

N ( bj_ 6f \ _ N < (6f - bj) 

"61-62 W - < &f - </ ~ (61 - 6a) - - <) ' 

For general n it will be useful to use the induction assumption (1791) also in the form 

n ~ 1 \N hN 



b N - b. 



S n (a (n) ,b (n) ) - S n (a (n) ) rr»-i. rr , l -2 / , v , v 



\ n-2 



with _, , . 

S B (a<»>) = -V"- 2 II^-O- ( 81 ) 

\/=l / ;<// 

In proceeding with the induction proof of (1501) we start separating terms involving 
the last summation (over p n -2) in (1791) and perform this sum using (1751) with f(t) = 
llm= 3 i (am -t): 



m 0"n—2 j 

1 ' rim=iK- fe m)y P t^ 2 "~" n"=i l K- 2 -6 s ) 



pi,...,p„_3 \i=l Llm=l\ 1 m > / p„_ 2 

( tt n?=f+i(^ -Qj) ^ 7y iv n m = i i( a m-6fc) 
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The first summation in the last equation involves also the phases of a m in the term 
n^li(Om ~~ fyfe) °f the last numerator. However, this term cancels just the term m = k 
in the denominator of the first bracket. So, comparing to ( 1791) for n — > n — 1 we can 
write 

SM n \ bW) = AT a- 2 £ - - &w -^- 2 B Li* } 7 , , (82) 



k= l \ a n-2 



where bj|. is the n— 2 component vector resulting from omitting the component 



Similarly, a^_ 2 is the n— 3-component vector resulting from a( n ) by omitting the 
component a n _2- Now we can use the induction assumption (180ft to insert here S* n _i: 

71—1 

S-nCaW.bW) = Na^_ 2 J2 ~^ :=r^ -x 



A:: 



[E ') n ngn?:W-<)' (83) 

It is more convenient to write (152"j) in terms of 5 n (a^ n ^) , equation (jSTl) . since then the 
b-dependence is explicit: 

"■— 1 iN I.N i 

k'<k" k k 11/=1 LU'=l(°l ~ a V ) 

X 



e n 



Apart from a sign, the terms linear in the by cancel since 

TTli- 1 / 7 7 \ ™ — l fe — 1 n— 1 

il?'^* _ k ' =k+1 k ' =1 s=l,s^k 

Using the same formula for the 6^-terms and simplifying the bf 1 — aff -terms, we get 



n-1 Tin-3( h N _ N\ Tl^fhN _ „iV 

Here the sum can be calculated using the identity 

n— 1 n— 3 n— 1 ,jy n— 3 

e nw-fln^ = -nw-«y. 

fe=i i=i *=i fc s 1=1 



n-2) 



s^k 

so that n - 3 



&(a<»>) = JVi a Vi(a£ a ) H( 



N „N 



a, —a 



n-2 



1=1 



which confirms (IHTj) and so also (IHU1) . □ 
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